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domain  model  of  the  pilot  and  a  least-squares  identification  algo¬ 
rithm.  NIPIP  functions  in  real-time  and  uses  a  ■"slidlng'^t  line  window 
to  maintain  freshness  in  the  data;  thus  time-varying  characteristics  in 
the  pilot's  control  strategy  can  be  measured. 

It  was  proposed  to  evaluate  this  technique  for  its  application  So  VTKS 
research.  STI  could  possibly  identify  pilot  behavioral  variations  as  a  !  unc¬ 
tion  of  task  changes  on  dependent  measures  of: 

1.  Pilot  input  bandwidth; 

2.  Pilot  stability  margin;  and 

3.  Crossfeed  control. 

In  particular,  development  of  proper  crossfeed  control  might  be  a  good  cri¬ 
terion  of  learning  for  glide  slope  control.  Hie  novice  pilot  is  unlikely  t ■>  he 
able  to  coordinate  power  and  pitch  adjustments  in  an  optimum  manner.  Hie  NIPIP 
program  may  be  able  to  identify  development  of  crossfeed  control,  or  any  break¬ 
down  in  the  strategy,  and  thus  could  provide  a  valuable  supplement  to  the 
existing  performance  measurement  package. 

The  first  set  of  data  supplied  to  STI  to  test  NIPIP  was  unsuitable  for 
complete  analysis  because  of  errors  in  the  turbulence  model  used  during  data 
collection.  More  data  were  collected  and  were  used  to  analyze  selected  runs 
from  an  aircraft  simulation  of  the  T-2C  on  final  approach  to  an  aircraft  car¬ 
rier.  The  NIPIP  results  presented  demonstrated  changes  in  the  pilot's 
describing  functions  with  simulated  glide  slope  disturbances  (injected  beam 
noise)  and  the  "tight"  versus  "loose"  tracking  runs.  For  the  "loose"  tracking 
runs,  there  was  a  very  low  glide  slope  gain  and  virtually  no  crossfeed  gain. 
For  the  "tight"  tracking  runs,  the  pilot  exhibited  high  glide  slope  and  cross¬ 
feed  gains  with  relatively  low  variability  in  the  data,  especially  for  tin*  runs 
with  beam  noise.  One  conclusion  is  that  the  NIPIP  technique  can  identify  the 
pilot's  control  strategy  and  behavior  by  parameters  meaningful  to  the  closed- 
loop  control  task.  Another  implication  is  that  adequate  glide  slope 
disturbances  must  be  present  in  order  for  the  pilot  to  demonstrate  his  ability 
to  control  the  aircraft  properly. 
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SECTION  I 
INTRODOCT ION 


This  report  contains  a  series  of  papers  prepared  by  Systems 
Technology,  Inc.,  (STI)  in  support  of  carrier-landing  research  in  the 
Visual  Technology  Research  Simulator  (VTRS).  The  following  work  was 
undertaken: 


1.  Development  of  a  quasi-random  turbulence  model.  This 
model  was  preferred  to  the  one  provided  initially  with 
the  VTRS  system  because  it  enabled  better  analysis  of 
pilot  responses  to  turbulence  inputs.  The  STI  model  is 
expected  to  be  appropriate  for  tasks  other  than  carrier 
landings  and  for  simulations  of  other  aircraft  types. 

2.  Modification  of  the  T-2C  simulation  to  more  closely 
represent  the  A-7  and  F-18  aircraft. 

3.  Application  and  evaluation  of  STI' s  Non-Intrusive  Pilot 
Identification  Program  (NIPIP),  which  was  developed  to 
estimate  the  pilot's  input-output  describing  function 
and  combined  pilot -vehicle  performance  parameters  such 
as  crossover  frequency  and  phase  margin  by  using  a  time 
domain  model  of  the  pilot  and  a  least-squares  identifi¬ 
cation  algorithm.  NIPIP  functions  in  real-time  and  uses 
a  "sliding"  time  window  to  maintain  freshness  in  the 
data;  thus  time-varying  characteristics  in  the  pilot's 
control  strategy  can  be  measured. 


It  was  proposed  to  evaluate  this  technique  for  its  application  to 
VTRS  research.  STI  could  possibly  identify  pilot  behavioral  variations  as 
a  function  of  task  changes  on  dependent  measures  of: 

1.  Pilot  input  bandwidth 

2.  Pilot  stability  margin 

3.  Crossfeed  control. 
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In  particular,  development  of  proper  crossfeed  control  might  be  a  good 
criterion  of  learning  for  glide  slope  control.  The  novice  pilot  is  un¬ 
likely  to  be  able  to  coordinate  power  and  pitch  adjustments  in  an  optimum 
manner.  The  NIPIP  program  may  be  able  to  identify  development  of  cross¬ 
feed  control,  or  any  breakdown  in  the  strategy,  and  thus  could  provide  a 
valuable  supplement  to  the  existing  performance  measurement  package. 

The  first  set  of  data  supplied  to  STI  to  test  NIPIP  was  unsuitable 
for  complete  analysis  because  of  errors  in  the  turbulence  model  used  dur¬ 
ing  data  collection.  More  data  were  collected  and  were  used  to  analyze 
selected  runs  from  an  aircraft  simulation  of  the  T-2C  on  final  approach  to 
an  aircraft  carrier.  The  NIPIP  results  presented  demonstrated  changes  in 
the  pilot's  describing  functions  with  simulated  glide  slope  disturbances 
(Injected  beam  noise)  and  the  "tight"  versus  "loose"  tracking  runs.  For 
the  "loose"  tracking  runs,  there  was  a  very  low  glide  slope  gain  and  vir¬ 
tually  no  crossfeed  gain.  For  the  "tight"  tracking  runs,  the  pilot 
exhibited  high  glide  slope  and  crossfeed  gains  with  relatively  low  vari¬ 
ability  in  the  data,  especially  for  the  runs  with  beam  noise.  The 
conclusions  are  that  a)  adequate  glide  slope  disturbances  must  be  present 
in  order  for  the  pilot  to  demonstrate  his  ability  to  control  the  aircraft 
properly,  and  b)  the  NIPIP  technique  can  objectively  and  quantitatively 
reveal  the  presence  or  absence  and  degree  of  a  specified  piloting  tech¬ 
nique  during  training  operations. 

Each  of  the  draft  reports,  submitted  in  working-paper  form  by  STI,  is 
given  In  a  separate  chapter  herein,  along  with  the  references  for  that 
chapter. 
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SECTION  II 

QUASI- RANDOM  COST  INPUTS  FOR  THE  T-2C  SIMULATION  AT  HTEC 


THE  PROBLEM 


The  training  scenario  is  simulated  carrier  landings  in  a  Buckeye  T-2C 
jet-trainer  onto  a  carrier  via  the  Fresnel  Lens  Optical  Landing  System 
(FLOLS).  A  number  of  simulator— device  variables  are  being  investigated, 
such  as:  raster-lines  (1025  versus  525),  brightness,  motion-base  (versus 
none),  haptic  seat  (versus  none),  sea-texture  (versus  none),  visual  scene 
lags  (approximately  0.10  sec  versus  0.20  sec),  pilot  experience  (basic- 
trained  versus  fleet-experienced),  several  FLOLS/carrier  representations 
(CGI  versus  model/CCTV,  each  with  fine  and  coarse  details),  and  others. 
The  carrier  deck  is  not  moving;  so  turbulence  inputs  are  used  to  provide  a 
relevant  control  task,  to  excite  the  closed-loop  vehicle  dynamic  modes 
under  pilot  control,  and  to  permit  later  analysis  of  selected  data  in  more 
detail. 

The  normal  simulator's  "carrier- burble"  is  not  used  during  this  ex¬ 
periment.  Hie  approach  takes  about  55  sec,  of  which  40  sec  are  on  the 
approach  path.  The  pilot  uses  throttle  to  control  airspeed  and  height, 
elevator  to  control  attitude,  and  aileron  to  control  lateral  lineup  er¬ 
rors,  according  to  standard  carrier  landing  procedures  (e.g.,  Refs.  1  and 
2*). 

The  gust  inputs  should  meet  the  following  criteria: 

1.  Quasi-random  appearing,  as  perceived  by  the  pilot,  i.e., 
subjectively  unpredictable  and  not  easily  or  unavoidably 
learned  upon  subsequent  encounters  (to  avoid  learning 
artifacts) . 


♦References  for  this  section  are  given  on  p.  11. 
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2.  Appreciable  frequency  content  near  the  frequency  regions 
involved  in  the  pilot-vehicle  loop  closures  ["phugoid," 
"path"  (height  and  offset),  and  "short-period"  and 
"Dutch-roll"  modes). 

3.  Repeatable  rms  levels  to  help  make  the  results  more 
consistent. 

4.  "Realistic"  and  justifiable  gust  spectra  for  the  turbu¬ 
lence  behind  a  carrier. 

5.  (Desirable)  - Permit  later  frequency  domain  analysis  to 

untangle  the  pilot's  control  strategy,  if  needed,  such 
as  crossfeeds  between  elevator  and  throttle,  effect  of 
motion,  etc. 

6.  Be  insertable  in  the  simulation  as  u  ,  v  ,  and  w  gusts. 

o  o  o 


THE  APPROACH 

A  set  of  three  independent  sum-of-sinusoid  inputs  is  ideal  for  the 
following  reasons: 


1.  By  keeping  integer  number  of  cycles  over  a  40  sec  run 

length,  each  rms  input  value  will  be  identical  from  run- 
to-run;  and  Fourier  analyses  will  show  distinct  line 
spectra;  while  the  random  initial  phase  will  provide 

different  appearing  waveforms. 

2.  If  all  frequencies  are  different  among  inputs,  then 

crossfeeds  from  one  axis  to  another  may  be  more  easily 
determined  (by  later  frequency-domain  analyses);  and  the 
inputs  are  statistically  independent. 

3.  The  frequencies  may  be  placed  to  excite  the  very  low 

frequency  closed-loop  modes  more  consistently  and  ef¬ 
ficiently  than  a  purely  random  signal. 

4.  The  amplitudes  may  be  shaped  to  match  the  effective 

power-spectral  density  of  actual  turbulence  in  the  de¬ 
sired  frequency  range. 
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A  check  of  Che  T-2C  dynamics,  based  on  data  given  in  Ref.  3  and 
elsewhere,  showed  Chat  Che  phugold  and  path  modes  lie  in  the  region  from 
«  ■  0.16  to  0.30  rad/sec  (0.03  to  0.06  Hz),  while  the  lateral  roll-spiral 
modes  lie  near  0.5  to  1.0  rad/sec,  the  pitch  short  period  near 
2.1  rad/sec,  and  Dutch-roll  mode  near  2  to  3  rad/sec  (depending  on  yaw 
damper  activation).  Therefore  inputs  spanning  0.03  to  3.0  rad/sec  are 
desired;  distributed  in  each  axis  as  appropriate  for  its  closed-loop  fre¬ 
quency  region.  This  requires  a  minimum  run  length  of  40  sec  for  which  the 
fundamental  frequency  is  0.025  Hz  or  0.157  rad/sec. 


RESULTS 

After  careful  consideration  of  non-simple  harmonic  ratios  in  each 
axis  (with  the  exception  of  the  1:3  radio  in  ug,  which  cannot  be  helped), 
integer  number  of  cycles  in  40  sec,  good  dispersion  on  the  log- frequency 
(Bode)  plots,  not-too-close  spacing  thereon,  and  no  more  than  4  to  5  sinu¬ 
soids  per  axis  for  good  signal/noise  considerations;  the  frequencies  shown 
in  Table  1  were  selected. 

• 

The  gust  spectral  shapes  were  based  on  the  actual  carrier  turbulence 
data  given  in  Ref.  2,  Appendix  B. 

For  this  special  case,  which  is  similar  to  but  not  quite  the  same  as 
low  altitude  atmospheric  turbulence,  the  spectral  density  shape  in  the 
vicinity  of  the  gust  input  region  was  closely  approximated  as  first-order 
filtered  noise.  To  compensate  for  the  non-uniform  spacing  of  sinusoids, 
each  sinusoid  is  assumed  to  cover  an  effective  bandwidth  set  by  the  geo¬ 
metric  mean  between  adjacent  frequencies,  and  an  equal  spacing  beyond  each 
last  one.  Given  the  input  frequencies,  phases,  and  desired  spectral  en¬ 
velope,  a  proprietary  computer  program  computes  the  required  amplitudes, 
amplitude-probability  distribution,  and  typical  waveform.  Samples  are 
Included  here  as  Figs.  1,  2,  and  3. 
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TABLE  1.  GUST  INPUT  SPECTRA  FOR  T-2C  LANDING  SIMULATION 


Number 
Cycles 
AO  sec 

of 

in 

Run 

Frequency 

N 

(Hz) 

“i 

(rad /sec 

1 

0.025 

0.157 

2 

0.050 

0.314 

3 

0.075 

0.471 

4 

0.100 

0.628 

5 

0.125 

0.785 

7 

0.175 

1.10 

8 

0.200 

1.26 

9 

0.225 

1.41 

11 

0.275 

1.73 

13 

0.325 

2.04 

14 

0.350 

2.20 

17 

0.425 

2.67 

19 

0.475 

2.83 

23 

0.575 

3.61 

Gust  Shaped  Amplitudes  for 


RMS  - 

1  unit  in  each 

Axis 

Comment  s 

Axial 

Lateral 

Vertical 

Modes 

(ft/sec) 

(ft/sec) 

(ft/sec) 

Excited : 

AUj  -  0.645 

aw1  -  i-175! 

.  Speed  and 
height  modes 

“K 

AUz  -  0.749 

AVi  =  1.295 

! 

\2  *  0.580  , 

lateral 

Path  Modes 

A,.  =  0.707 
^*3 

1 

AW3  -  0.381  | 

|i/t’2 

Au4  * 

A„  »  0.431 

2 

} 

1 

\  -  °-254  | 

Short 

Period  Mode 

to 

1  sp 

=■  0.288 
^3 

i 

V.  "  0-473 

Av  -  0.234 

V4 

1 

AW5  -  0.226  | 

i 

Dutch 

Roll 

Mode 

,JDR 

RMS  -  I  A^/2 
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Each  signal  spectra  has  been  adjusted  for  1  unit  (e.g.,  ft/sec)  rms 
value,  and  the  sun  should  be  multiplied  by  the  desired  turbulence  inten¬ 
sity  levels.  For  this  carrier  landing  case,  use  the  same  rms  gust 
intensity  in  each  axis,  as  follows: 


GUstiness 

RMS  Level 
f  t/sec 

Low 

1.0 

Moderate 

3.0 

Rough 

5.0 

We  recommend  3.0  ft/sec  as  a  representative  level.  Note  that  peak  gust 
velocities  will  occasionally  reach  3  to  10  ft/sec.  This  will  not  stall 
the  wing  or  tails,  yet  is  big  enough  to  provide  adequate  path  and  motion 
system  excitation. 

Rather  than  provide  purely  randomized  mutual  phasing,  we  recommend 
using  3  or  4  sets  (of  phases),  so  that  corresponding  time  traces  may  be 
overlaid  for  comparisons.  Each  of  the  sets  should  be  checked  for  a  sub¬ 
jectively  "random"  appearance,  i.e.,  no  large  waveform  "signatures"  (which 
are  detectable  upon  third  encounter) .  They  can  be  used  in  random  or  sys¬ 
tematic  order  to  avoid  learning  artifacts. 

Even  though  the  range  to  the  FLOLS  is  varying  during  the  approach, 
experience  has  shown  that  fairly  stationary  statistics  exist  along  the 
approach  path  until  the  last  5  to  10  sec.  Therefore  the  approach  perfor¬ 
mance  rms  statistics  should  be  based  on  the  time  period  from  50  to  10  sec 
before  nominal  touchdown,  and  separate  (transient-type)  statistics  should 
be  used  for  the  last  10  sec  (as  well  as  for  the  first  5  sec  after  the 
start  of  final  approach). 
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SECTION  III 

REVIEW  OV  THE  NTEC  T-2C  LANDING  SIMULATOR  MATHEMATICAL  MODEL 
AMD  ITS  MODIFICATION  TO  A  FLEET-LIKE  FIGHTER 


INTRODUCTION 

Simulator  tests  of  a  modified  FLOLS  visual  landing  aid  (the  rate-error 
display  bars  of  Lt.  C.  Kaul)  are  currently  being  planned  for  the  simulator 
at  the  Naval  Training  Equipment  Center  in  Orlando,  Florida.  The  aircraft 
mathematical  model  was  reviewed  with  a  view  to  modifying  its  characteris¬ 
tics  to  resemble  current  Navy  fighters,  thereby  providing  a  carrier 
approach  task  more  representative  of  fleet  operations  which  might  benefit 
from  the  improved  FLOLS  display.  The  major  data  source  is  an  NTEC  FORTRAN 
listing  of  the  "Aerodynamic  Subroutine"  which  models  the  T-2C,  supple¬ 
mented  by  an  antecedent  dociment,  the  Singer  Link  T-2C  simulation 
subroutines  (Ref.  1*). 

MODEL  REVIEW 

The  NTEC  FORTRAN  listing,  dated  17  July  i980,  was  reviewed  and  com¬ 
pared  against  the  Sinker  Link  routines  in  Ref.  1;  and  the  resulting  "non- 
dimensional"  aerodynamic  force  and  moment  equations  were  abstracted  in  the 
format  of  Ref.  2.  The  following  anomalies  were  noted: 


1. 


The  sign  of  Z$  (due  to  Cl,  )  appears  incorrect.  For 
e  Og 

aft  horizontal  tail,  deflection  of  the  elevator  (e.g. , 
Tg  up)  should  result  in  the  tail's  pitching  moment 


♦References  for  Section  III  are  on  p.  26. 
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(positive  nose  up)  and  the  z-axis  force  (positive  down¬ 
ward)  being  of  the  saute  sign. 

2.  The  sign  of  Cj^  appears  incorrect.  For  subsonic  air¬ 
craft,  the  sign  of  this  term  should  be  the  same  as  . 
Also,  a  modification  has  been  Inserted  in  the  code  mul¬ 
tiplying  this  coefficient  by  0.47. 

3.  The  combination  of  +  0.47C^a]  is  multiplied  by  a 

factor  of  4.0,  resulting  in  a  computed  value  for 
approximately  4  times  that  given  for  the  T-2C  In 
Table  11  of  Ref.  3  for  a  similar  flight  condition. 


It  is  hypothesized  that  these  errors  tend  to  cancel  each  other  from 
the  viewpoint  of  the  pilot  in  the  cab.  Certainly  increased  aircraft  damp¬ 
ing  would  tend  to  mitigate  the  otherwise  abrupt  normal  acceleration 
response  to  elevator  deflection.  It  would  also  tend  to  negate  the  desta¬ 
bilizing  tendency  associated  with  the  positive  sign  on  C^* 

The  errors  were  corrected  and  dimensional  coefficients  were  calculated 
for  comparison  with  the  data  of  Ref.  3  and  the  wind  tunnel  data  cited  in 
Ref.  4.  This  comparison  is  shown  in  Table  1.  Experience  relating  various 
derivatives  to  response  dynamics  (e.g..  Ref.  2)  suggests  substantial 
agreement  for  the  derivatives  affecting  high-  and  mid-frequency  re¬ 
sponses.  Some  differences  are  apparent  in  the  derivatives  affecting  * ow- 
frequency  speed-related  responses,  which  may  be  attributable  to 
differences  in  flaps,  speed  brake,  or  landing  gear  positions. 
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TABLE  1.  T-2C  POWER  APPROACH  DEVIATION  COMPARISON 


Parameter 

Units 

"Corrected" 

Simulation 

Ref.  3 

Ref.  4 

Velocity 

(kt) 

ft/sec 

(108) 

183 

(108) 

183 

(139.6) 

236 

Weight 

lb 

11,000 

11,000 

1 ,000 

“trim 

deg 

6.99 

.  * 

4.0 

4.7 

Flap 

deg 

33. 

— 

16. 

Speed 

brake 

deg 

0 

— 

•  0 

Gear 

Down 

— 

Down 

sec”1 

-  0.0742 

-  0.046 

-  0.046 

Zu 

sec-1 

-  0.3427 

-  0.2206 

-  0.269 

Mu 

(ft-sec)”1 

0.00092 

— 

0.00434 

Mw 

ft”1 

-  0.002249 

-  0.002 

-  0.00225 

X« 

sec”1 

-  0.1350 

0.0534 

-  0.0983 

*w 

sec”1 

-  0.7384 

-  0.7683 

-  0.974 

“w 

(f t-sec)”1 

-  0.01531 

-  0.0168 

-  0.01945 

sec”1 

-  1.1593 

-  1.105 

-  1.42 

X«e 

ft/secZ 

-  6.987 

— 

— 

Z«e 

f t/secZ 

-15.54 

014.45 

-24.07 

M«e 

sec”Z 

-  5.187 

-  5.824 

-  9.63 

x«T 

f  t/lb-secZ 

0.002904 

— 

— 

f t/lb-secZ 

-  0.000352 

— 

— 

M«t 

(lb-secz)  1 

0.0000363 

— 

— 

*Based  on  stated  value 
stability  axes  of  +0.07  rad 

of  thrust  line  inclination  relative 
(+4.0  deg). 

to 
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liable  2  lists  some  of  the  sore  important  transfer  function  parameters 
associated  with  the  corrected  simulation  data*.  Comparison  with  handling 
qualities  experience  (e.g. ,  Kef.  2)  shows  the  airplane  to  be  "well  be¬ 
haved"  in  its  attitude  responses  (well-damped  short  period  near 
2  rad/sec),  slightly  on  the  backside  (small  negative  value 
for  l/T^j),  with  adequate  path  angle  response  bandwidth  ( 1 / Tq2  >  0.7). 
Because  the  throttle  response  lag  for  the  T-2C  is  relatively  short,  on  the 
order  of  0.25  sec,  the  path  response  to  the  throttle  is  good,  as  is  the 
path  response  to  attitude  change.  The  aircraft  also  exhibits  a  modest 
pitch  up  with  power  application. 

Either  the  conventional  CTOL  (primary  emphasis  on  attitude  change  for 
path  correction)  or  the  "Navy  doctrinal"  (primary  emphasis  on  throttle  for 
path  correction)  techniques  are  usable  for  this  airplane,  the  former  with 

some  modification  to  correct  for  "backsidedness"  -  being  on  the  "back 

side"  of  the  power- required  curve  (i.e.  ,  below  minimum-power  speed). 


MODEL  MODIFICATIONS 

To  make  the  T-2C  behave  more  like  a  typical  Navy  fighter  on  carrier 
approach  or,  more  particularly,  to  make  it  behave  like  an  aircraft  known 
to  have  less  than  desirable  handling  qualities  for  this  task  requires  that 
the  path  responses  to  throttle  and  elevator  be  determined.  The  A-7 
(Ref.  5)  and  F-18-like  (Ref.  6)  aircraft  provide  good  examples.  The  chan¬ 
ges  (relative  to  the  "corrected  simulation"  of  the  T-2C)  required  to 
deteriorate  the  responses  are  as  follows: 


*A  common  shorthand  format  is  used  to  list  the  first-  and  second-order 
dynamics: 

high-frequency  gain  =  xxx;  (s  +  1/T)  =  (1/T);  [s2  +  2Cws  +  id21  =  U,u>l 
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TABLE  2.  T-2C  TRANSFER  FUNCTION  FACTORS  FOR 
POWER  APPROACH  (CORRECTED  DERIVATIVES) 

Format  ia:  Gain  (I/T) 


Denominator 


A  -  [0.0653,  0.200] [0.614,  1.919] 

Cap  “^p  C8p 

Elevator  Numeratora 


0 

Nsfi  -  “5.15(0.01552X0.756) 

1/Tflj  1/Tq2 

N5e  *  15. 54( -0.0695)  (—5. 77)  (7 .33) 

1/  Td  j  l/Tdj 

N6e  "  “6.99(0. 420)  (-5. 54)  (7. 13) 
l/TUj  1/TU2  1/TUj 

Throttle  Numeratora 


N6t  -  0.000352(0.303) [0.466,  4.48] 
N«T  -  0.00290(-0. 1106) [0.678,  1.796] 

Coupling  Numerators  (Throttle  and  Elevanr 

N6T6e  "  -°*01470(0.716) 

A  0 


0.00239(2.18) 


» 
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1.  Airspeed  of  130  kt  (219  ft/sec)  -  like  the  A-7. 

2.  Flaps  to  16  deg,  speed  brakes  fully  extended,  gear  down. 

3.  Lift  equation  changes  (provides  required  Cg  and  3C^/3a 
at  the  desired  value  of  a): 

AFKCLALP  -  0.65 
AFKCLFLP  -  1.38 
AFKCLWKR  -  0 

4.  Drag  equation  changes  (provides  required  CD  and  SCp/Sa): 

AFKCDCL  -  1.51 
AFKCDSB  -  1.18 
AFKCDFW  -  0 

5.  Engine  equivalent  time  constant,  Tg  ■  1.0  sec. 

Table  3  lists  the  resulting  stability  derivatives  and  compares  them  with 
the  A-7  and  F-18-like  data  of  Refs.  5  and  6.  Appended  to  the  table  is  a 
selected  set  of  transfer  function  factors.  These  data  suggest  that  the 
path  responses  would  be  similar  for  the  three  cases.  The  short-period 
response  differences  are  not  so  great  as  indicated  because  both  the  A-7 
and  F- 18-1 ike  aircraft  are  modified  by  their  flight  control  systems  to 
revise  their  effective  short— period  frequencies. 

With  regard  to  its  path  responses  to  either  the  stick  ^att'!u-^c 
change)  or  the  throttle,  the  modified  T-2C  should  exhibit  properties  ‘,:mi- 
lar  to  those  of  the  A-7  and  F-18.  In  brief,  the  pilot  will  find  It 
difficult  to  fly  on  approach  because  its  responses  simply  are  not  fast 
enough  [slow  engine,  long  (2+  sec)  path  response  lag,  T@2].  Further,  the 
closure  rate  is  similar  to  that  of  the  other  two  aircraft,  which  should 
present  the  pilot  with  similar  time-to-go  stress  and  perceptual- range 
difficulties  (Ref.  7). 

All  three  aircraft  suffer  from  deficient  responses  in  the  "outer" 

control  loop  for  precision  path  control  in  carrier  approach  _  deficient 

in  the  sense  that  they  are  not  fast  enough  and  (possibly  barring  the  F-18 
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TABLE  3.  DERIVATIVE  COMPARISON,  MODIFIED  T-2C  SIMULATION 


WITH  FLEET  AIRCRAFT 

Parameter 

Units 

Modified  T-2C 

A-7 

F-18-Like 

Simulation 

(Ref.  5) 

(Ref.  6) 

Velocity 

ft/sec 

219 

218 

226 

Weight 

lb 

13,200 

24  ,000 

27  ,890 

“trim 

deg 

6.9 

12.0 

5.5 

see"* 

-  0.0608 

-  0.0545 

-  0.0691 

sec-1 

-  0.2865 

-  0.2870 

-  0.2950 

«u 

(ft-sec)'1 

-  0.000755 

-  0.000165 

0 

ft-1 

-  0.002249 

-  0.000289 

-  0.0003167 

sec'1 

0.0674 

0.0643 

0.0777 

Zw 

sec'1 

-  0.6008 

-  0.5289 

-  0.4822 

Mw 

(f t-sec)-1 

-  0.01221 

-  0.007964 

-  0.001677 

Mq 

sec'1 

-  1.387 

-  0.3275 

-  0.2333 

X«e 

(f t-sec) Z 

-  3.379 

0.7328 

4.549 

Z«e 

f  t/sec^ 

-22.16 

-14.71 

-17.43 

M«e 

sec'Z 

-  7.141 

-  2.189 

-  1.75 

f t/lb-secZ 

0.002903 

0.001317 

0.00148 

Z 6 x 

f t/lb-secZ 

0.000353 

-  0.000250 

-  0.0001106 

Mfix 

,  -1 

(lb-sec^) 

0.000036 

0.000004 

0. 

a) 

ap 

rad /sec 

1.88 

1.38 

0.71 

^sp 

i/t02 

1/sec 

0.522 

0.429 

0.397 

1/Tdl 

1/sec 

0.0232 

0.0030 

0.0246 

*E 

sec 

1.0 

1.0 

0.54 
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with  its  faster  engine)  that  they  force  the  pilot  to  make  path  corrections 
with  the  elevator  when  throttle  response  delay  becomes  critical. 

All  three  aircraft  will  benefit  if  low-frequency  path  deviation  rate 
("lead  equalization")  Is  introduced  in  the  outer  loop.  The  modified  FLOLS 
scheme  offers  this  potential  by  providing  the  pilot  with  an  explicit  Indi¬ 
cation  of  path  error  rate.  We  would  expect  to  see  improved  pilot  opinion 
(due  to  workload  reduction)  and  perhaps  even  Improved  performance  with 
such  a  scheme.  The  improvement  should  be  more  noticeable  with  more  dif¬ 
ficult  aircraft. 


PILOT  LOOP  CLOSURE  TOPOLOGIES 

Among  the  measurements  which  have  been  proposed  for  this  experiment  is 
the  determination  of  pilot  describing  functions  In  the  longitudinal  ap¬ 
proach  task.  These  will  endeavor  to  establish  pilot  gain  and  equalization 
in  the  several  loops,  hence  the  control  technique  being  employed.  The 
material  which  follows  describes  the  loop  closure  possibilities. 

Common  to  all  is  control  of  pitch  attitude  with  elevator,  the  percep¬ 
tual  cue  coming  from  the  pilot' s . awareness  of  where  the  aircraft  nose  is 
relative  to  the  horizon,  or,  failing  this  (e.g.,  at  night),  from  the  car¬ 
rier  lights  themselves.  If  this  loop  is  closed  loosely,  subsequent 
closures  of  path  deviation  to  either  throttle  or  attitude  are  limited  to 
low  frequencies  in  the  vicinity  of  the  phugoid  mode.  If,  on  the  other 
hand,  the  attitude  loop  is  closed  tightly,  it  becomes  possible  to  achieve 
higher  path-control  response  bandwidths,  limited  by  l/Tg^  *  -Zw.  This  can 
be  accomplished  with  outer-loop  path  deviation  closures  to  either  the 
attitude  command  or  to  the  throttle,  provided  that  the  pilot's  control 
technique  corrects  for  the  airspeed  deviations  otherwise  incurred.  Since 
he  does  not  have  a  head-up  indication  of  speed  (both  his  angle-of-attack 
indexer  and  the  perceived  closure  rate  to  the  carrier  are  poor  indicators 
of  airspeed),  the  skilled  pilot  generally  accomplishes  this  via  "cross¬ 
feeds"  to  the  other  control. 
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Consider  first  the  preferred  ("Navy  doctrine")  technique  -  path  con¬ 

trol  via  the  throttle  with  a  tight  attitude  closure  via  elevator. 
Figure  la  shows  the  topology  of  the  control  loop  structure  with  the 
addition  of  a  crossfeed  from  the  throttle  to  the  attitude  command.  With  a 
tight  attitude  loop  closure  the  phugoid  mode  is  driven  into  the  attitude 
numerator  zeros,  the  short-period  mode  to  higher  frequencies  such  that  the 
path  deviation  response  to  thrust  changes  can  be  approximated  by  (thrust 
lag  ignored) : 


*d  9 
6_6 
T  e 


(s)N®  [C'  u »'  ] 
O  sp  sp 
e 


which,  for  the  T-2C,  is: 


_ -0.00239(2.18) _ 

-5. 15(0)(0. 01552)  (0.756)  [  c'  ,w'  ] 

sp  sp 

1/Tg  1/ Tg  Closed-pitch 

1  2  loop  dynamics 


Clearly  it  turns  out  (Ref.  5)  that  the  ability  of  the  pilot  to  control 
path  is  limited  by  the  lag  at  1/Tq^  or,  more  precisely,  by  the  closed-loop 
root  resulting  from  0  ♦  6e. 

Now  suppose  that  the  pilot  is  skilled  enough  to  pull  the  stick  so  as 
to  raise  the  nose  slightly  upon  power  application  (this  tends  to  occur 
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a.  Navy  Technique,  Crossfeed  for  Speed  Control 


b.  Conventional  Technique,  Crossfeed  for  Speed  Technique 


Figure  1 .  Loop  Structures  for  Carrier  Approach 
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naturally  for  the  T-2C  with  Its  pltch-wlth-power  characteristic)  such  as 
to  minimize  the  resulting  speed  perturbation.  That  is: 


9 


n“  «9 

5_  e 

_ I _ +  Y  e 

xi0  5 

N,  e 


0 


This  requires  a  crossfeed  of: 


^  6 
T  e 


n: 


•  -0.01470(0.716) 

-6. 99(0 .420)  (-5. 54)  (7 . 1 3) 


which,  at  frequencies  below  0.42  rad/sec,  Is  equivalent  to  a  gain 
0 

Y^c  ■  0.00009076  rad  (0c)/lb  (5^).  The  resulting  path  response  is  given 
by: 


d_ 

6. 


5  6 
T  e 


9 


9+5 

6_-*-8e 
T  c 


(s)N.  [?',<) 
o  sp  sp 


‘5t  (s)N®  U'  ] 
sp  sp 


The  effective  numerator  is: 


d  9 

V  + 

T  e 


9  d 
Y6Cn5 


9 

-0.00239(2.18)  +  Y.C  15.  54(-0.0695)  (-5.  77)  (7 .33) 
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which  can  be  approximated  at  low  frequencies  by: 


Ns  «  + 

T  e  T  e 


-0.06204(0.01716) 


The  resulting  response  is: 


Approximately 
cancels  1/T« 

61 

-0.06204(0.01716) 

-5.  15(0)  (0.01552)  (0.7  56)  [c'  ,<*>'  ] 

sp  sp 


1/T  1/T0 

1  2 

Now  the  dominant  performance-limiting  lag  is  at  1/Tq2,  considerably  higher 
than  before. 

Some  remarks  and  caveats: 


1.  Cutting  the  crosafeed  gain  to  about  60  percent  of  the 

value  given  here  shifts  the  zero  to  the  vicinity  of 
0.07,  which  represents  a  more  realistic  value 

of  1/T^j,  the  closed-loop  speed  response  dynamics. 

2.  Thrust  response  lag  will  limit  the  bandwidth  of  the 
d  ♦  6m  closure.  Here  it  has  been  improved  because 
Tg  ■  0.25  sec  for  the  T-2C,  considerably  faster  than  the 
2  sec  (or  thereabouts)  closed-loop  path  response  lag 
obtainable  with  this  technique. 

3.  While  speed  perturbations  are  minimized,  it  will  still 
be  necessary  to  close  a  u  ♦  loop  intermittently  to 
counter  the  effects  of  the  carrier  wake  "burble.”  Ac¬ 
tually,  skilled  pilots  may  even  act,  to  a  degree, 
"precognitively"  on  this;  i.e.,  knowing  the  carrier  and 
wind  over  deck  speed  and  direction  and  then  range,  they 
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put  in  anticipatory  throttle  and  elevator  control  pulses 
to  compensate  for  the  burble  as  they  pass  through  It. 


Attention  is  now  directed  to  the  more  conventional  CTOL  technique 

represented  by  the  loop  structure  of  Fig.  lb.  Here  again  a  crossfeed  is 

introduced,  this  tine  to  the  throttle.  It  has  the  saae  purpose  and,  in 

0 

fact,  is  given  by  the  inverse  of  the  of  the  previous  example.  The 
crossfeed  moves  l/T^,  the  low-frequency  (right  half  plane  in  this  case) 
zero  of  the  N$T  nuaerator  to  the  left  until  it  "cancels"  1/ T§ ^ ,  thereby 
allowing  d  ♦  8C  closures  of  bandwldths  approaching  l/T^,  limited  in  this 
case  by  the  short-period  response  dynamics  of  the  aircraft. 

In  examining  Fig.  1  in  the  context  of  the  forthcoming  experiment  we 
note  the  following: 


1.  Describing  function  measurements  of  pilot  behavior  will 
reveal  the  existence  of  path  deviation  error-correlated 
signals  on  both  throttle  and  elevator  _if  the  speed  re¬ 
sponse  cancelling  crossfeed  is  in  operation  with  either 
loop  structure.  Only  if  the  crossfeed  is  absent  can 
such  measurements  distinguish  between  the  two  techniques 
for  path  control. 

2.  Flight  records,  e.g.  ,  Fig.  7  of  Ref.  5,  suggest  time- 
varying  pilot  behavior  characterized  by  increasing  pilot 
gains  as  the  carrier  ramp  is  approached.  Thus  it  is 
quite  likely  that  simpler  loop  structures,  having  less 
performance  potential,  are  used  early  in  the  approach  to 
minimize  workload,  while  the  higher-performance  loop 
structures  and  gain  are  used  later  when  needed.  The 
time-varying  describing  functions  are  anticipated  to 
show  time-varying  behavior  as  the  ramp  is  approached. 

3.  The  availability  of  outer-loop  lead  compensation  from 
the  modified  FLOLS  should  ease  the  pilot's  task  with 
either  loop  topology,  particularly  when  the  path  re¬ 
sponse  is  limited  by  one  or  more  of: 

a.  Low  1/Tg2 


b.  Slow-responding  engine 
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Thus  we  would  anticipate  that  the  "modified  T-2C"  of 
Table  3  might  show  greater  improvement  from  rate  bars  in 
pilot  opinion,  and  perhaps  even  performance,  than  the 
T-2C,  per  se. 

4.  Speed  deviations  caused  by  external  gusts  and  turbulence 
are  anticipated  to  be  mostly  correlated,  if  at  all,  with 
the  throttle  and  stick  during  the  final  seconds  of  the 
approach.  The  pilot  cannot  perceive  such  errors 
directly;  rather,  he  can  only  respond  to  the  resulting 
path  errors  and  to  the  changing  time  indication  provided 
by  his  indexer.  A  disturbance  artificially  introduced 
on  the  Indexer  alone  could  reveal  closed-loop  usage  of 
this  Indicator. 


This  analysis  has  pointed  out: 


1.  Some  minor  discrepancies  in  the  T-2C  simulator  mathe¬ 
matical  model  which  should  be  resolved  before  this 
simulation. 

2.  Simulator  modifications  which  will  result  in  path  re¬ 
sponses  (not  short-period  responses)  more  typical  of 
operational  Navy  fighters. 

3.  Anticipated  signal  correlations  and  trends  in  pilot 
describing  function  measurements. 


25 


NAVTRAEQUIPCEN  78-C-0060-10 


REFERENCES  FOR  SECTIOH  III 


1.  Anon.,  Visual  Technology  Flight  Simulator  .VxrS)  Math  Model  Report, 

Singer-Link  Division,  NAVTRAEQUIPCEN7 5-C-0009-1 1,  March  1977. 

2.  McRuer,  Duane,  Irving  Ashkenas,  and  Dunstan  Graham,  Aircraft  Dynamics 

and  Automatic  Control,  Princeton  University  Press,  1973. 

3.  Smith,  R.  M. ,  and  N.  D.  Geddes,  Handling  Quality  Requirements  for 

Advanced  Aircraft  Design:  Longitudinal  Mode,  AFFDL-TR-78-154, 

August  1979. 

4.  Buenz,  D.  A.,  et  al. ,  Identification  of  T-2  Aerodynamic  Derivatives 

from  Flight  Data,  NADC  Report,  March  1975  (AD-A021116) . 

5.  Craig,  S.  J. ,  R.  P.  Ring land,  and  I.  L.  Ashkenas,  An  Analysis  of  Navy 

Approach  Power  Compensator  Problems  and  Requirements,  Systems 
Technology,  Inc.,  Technical  Report  No.  197-1,  March  1971. 

6.  Rlngland,  R.  F.,  and  D.  E.  Johnston,  Analytical  Assessment  of  the 

F-18A  Flying  Qualities  During  Carrier  Approach,  Systems  Tech¬ 
nology,  Inc.,  Technical  Report  No.  1090-1,  Sept.  1977. 

7.  Clement,  W.  F.,  and  R.  K.  Heffley,  An  Example  of  a  Concept  for  Es¬ 

tablishing  Flight  Training  Media  Requirements.  Systems  Technology, 
Inc.,  Technical  Report  No.  2108-1,  May  1980. 


26 


NAVTRAEQUIPCEN  78-C-0060-10 


SECTION  IV 

APPLICATION  OP  THE  HON- INTRUSIVE  PILOT  IDENTIFICATION 
PROGRAM  TO  A  MULTI- LOOP  CONTROL  TASK 

IHTRODOCTION  AND  BACKGROUND 

The  Navy  carrier  approach  control  task  is  analyzed  using  the  Non- 
Intrusive  Pilot  Identification  Program  (NIPIP,  Refs.  1,  2,  and  3*).  The 
result  is  the  identification  of  the  basic  pilot  control  strategy  in  con¬ 
ventional  feedback  control  terms. 

The  final  approach  to  an  aircraft  carrier  using  the  Fresnel  Lens 
Optical  Landing  System  (FLOLS)  is  analyzed  in  terms  of  two  loops:  flight 
path  (FLOLS)  control  and  speed  (angle  of  attack)  control.  The  nominal 
control  strategy  used  by  the  pilot  is  depicted  in  the  block  diagram  of 
Fig.  1.  The  objective  of  the  pilot  is  to  regulate  FLOLS  glide  slope  de¬ 
viation,  e,  and  angle  of  attack,  a,  against  external  disturbances  due  to 
both  axial  gusts,  u  ,  and  vertical  gusts,  w  .  The  e  and  a  loops  are  com- 

©  O 

monly  referred  to  as  "outer  loops."  The  pilot  must  also  regulate  pitch 
attitude  with  the  elevator,  which  is  an  "inner  loop." 

The  prescribed  Navy  piloting  technique  for  controlling  the  aircraft 
Is  to  regulate  e  with  the  throttle,  6^,  and  a  with  the  elevator,  5g  (via 

changes  in  attitude,  9) - a  so-called  "backside"  piloting  technique.  In 

reality,  however,  a  pilot  learns  that  he  must  also  "crossfeed"  the  throt¬ 
tle  to  pitch  attitude  in  order  to  achieve  adequate  response.  That  is, 
when  the  pilot  makes  a  correction  to  e  using  5^  he  also  adjusts  0X,  as 
shown  in  Fig.  1.  With  practice  a  pilot  will  learn  how  much  crossfeed  to 
use  for  a  given  aircraft  and  approach  flight  condition. 


♦References  for  Section  IV  are  on  p.  55. 
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Reference  4  demonstrates  that  the  pursuit-crossfeed  piloting  tech¬ 
nique  described  above  must  be  used  in  order  to  obtain  adequate  pilot- 
vehicle  performance.  Furthermore,  once  the  pilot  is  established  at  the 
reference  angle  of  attack,  the  crossfeed  significantly  lessens  the  need 
for  the  a  >  6g  loop  shown  in  Fig.  1.  The  purpose  of  this  report  is  to 
show  hew  the  control  technique  depicted  in  Fig.  1  can  be  quantified  using 
NIPIP.  Specifically,  the  objective  is  to  quantify  the  elements  of  the 
pilot's  control  technique,  i.e.  ,  Yp,  y£,  and  Yp . 

NIPIP  quantifies  the  terms  Yp,  Yp,  and  Yp  in  terms  of  frequency- 
domain  describing  functions  (Ref.  5).  The  pilot  describing  functions 
reflect  pilot-behavioral  performance  and  can  be  used  to  assess  pilot  work¬ 
load.  The  pilot-alone  performance  is  sometimes  a  more  sensitive  measure 
than  the  combined  pilot-vehicle  performance  (e.g.,  rms  glide  slope  devia¬ 
tion,  rms  airspeed  error,  etc.)  because  the  pilot  adapts  to  changes  in 
experimental  conditions  (e.g.,  vehicle  dynamics,  display  symbology,  etc.) 
in  order  to  maintain  the  combined-vehicle  performance  at  some  acceptable 
level.  For  example,  the  gradual  learning  of  the  proper  "crossfeed"  be¬ 
havior  between  the  throttle  and  the  elevator  inputs  can  be  measured  as  an 
index  of  skill  acquisition. 


DESCRIPTION  OF  THE  EXPERIMENT 

A  piloted  aircraft  simulation  was  performed  on  the  Visual  Technology 
Research  Simulator  (VTRS)  at  the  Naval  Training  Equipment  Center  (NTEC)  at 
Orlando,  Florida.  The  experiment  was  conducted  by  personnel  from  the 
Canyon  Research  Group,  Inc.,  using  an  experienced  Navy  pilot.  The  experi¬ 
mental  scenario  was  as  foll<x*s: 

1.  Night  approach  to  an  aircraft  carrier  using  a  computei~- 
generated  image  (CGI)  display 

2.  Raw  FLOLS  display 
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3.  Constant  rms  turbulence  using  the  sum  of  sine  waves 

technique  described  in  Table  1.  Only  the  vertical  (“g) 
and  longitudinal  (u  )  components  were  used  in  this  ex¬ 
periment  because  the  lateral  component  (v  )  was  used  to 
simulate  beam  noise.  The  rms  level  for  Both  u„  and  w 
was  3.0  fps.  g 

4.  No  aircraft  carrier  motion;  the  beam  noise  provided  a 
surrogate  for  the  carrier  motion  artifacts. 


The  experimental  matrix  consisted  of  variations  in  the  following 
parameters: 


1.  Constant  and  variable  range  from  aircraft  carrier  (l.e., 
constant  range  means  that  the  aircraft  was  not  allowed 
to  move  closer  to  the  aircraft  carrier)  .  Three  constant 
range  and  one  variable  range  conditions  were  used: 

a.  Variable  (normal  approach  starting  from  an  initial 
range  of  9000  ft) 

b.  Constant  "long"  range  (7200  ft) 

c.  Constant  "medium"  range  (3600  ft) 

d.  Constant  "short"  range  (1800  ft) 

2.  With  and  without  the  simulated  beam  noise  defined  in 
Table  2. 

3.  "Loose"  and  "tight"  tracking  of  the  glide  slope 
(Student  0  versus  Student  2  in  Table  3). 


The  actual  experimental  matrix  is  defined  in  Table  3.  The  experi¬ 
mental  conditions  were  not  randomized. 

The  T-2C  aircraft  was  used  for  all  of  the  runs  shown  in  Table  3,  and 
all  runs  were  flown  by  the  same  pilot.  However,  on  the  "Student  0"  runs, 
the  pilot  was  allowed  to  fly  the  (simulated)  aircraft  as  he  saw  fit.  The 
result  was  that  he  did  not  actively  track  the  glide  slope,  and  admitted 
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TABLE  1.  GUST  INPUT  SPECTRA  POR  T-2C  LANDING  SIWLATION* 


Nun  bar 

of 

Cycles 

la 

Gusc  Shaped  Anplltudes  for 

40  sec 

Run 

Frequency 

RMS  - 

lurvlt  In  each  Axis 

Comment s 

N 

“l 

Axial 

Lateral  t 

Vertical 

(rad /sec) 

(ft/sec) 

(ft/sec) 

(f t/aec) 

1 

0.025 

0.157 

*u,  - 

Speed  and 

1 

height  modes 

2 

0.050 

0.314 

\  -  l'175 

% 

3 

0.075 

0.471 

*u,  -  °-749 

2 

Letera  1 

4 

0.100 

0.628 

AVi  -  1.295 

Path  >todes 

S 

0.125 

0.785 

A_  -  0.580 
'*2 

7 

0.175 

1.10 

v,  •  °-707 

l/?e 

8 

0.200 

1.26 

Av^  -  0.381 

z 

9 

0.225 

1.41 

Av  -  0.431 

2 

Short 

11 

0.275 

1.73 

\  -  °-548 

Period  Made 

13 

0.325 

2.04 

AV(,  -  0.294 

Ui 

»p 

14 

0.350 

2.20 

A^  -  0.288 

17 

0.425 

2.67 

Au  -  0.473 

Dutch 

J 

Roll 

19 

0.475 

2.83 

Aw  -  0.226 

5 

Node 

23 

0.575 

3.61 

-  0.234 

4 

sa|  -  2.00 

2.0 

2.0 

RMS  -  (lA7/2 

)1/2  -  1.0  for 

“8-  V  “g 

Each  co<a>onent  of  turbulence  la 

calculated  at  follows:  ■ 

I,  1  *1  Sin  (uijt 

*  *1 

where 

**  • 

V  V  or 

V 

K*  •  seal*  (actor  •  3.0  for  thia  slaulatlond  .e.  m  gust  level  la  3.0  (pa); 

and  are  randoa  phaaa  anglea  that  are  constant  throughout  each  run  but  change 

froa  run  to  run. 


‘Adapted  froa  Raf.  6. 

7For  this  experlaenc  only  the  lateral  coi^onents  are  deleted  and  are  used  for  bean 
noise  (see  Table  2). 
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TABLE  2.  FLOLS  BEAM  NOISE 


Nunber  of  Cycles 

In  40  sec  Run  Frequency  Amplitudes 


“i. 

Ae< 

(Hz) 

(rad/sec) 

i 

(rad) 

4 

0.100 

0.628 

0.003503 

9 

0.225 

1.41 

0.003702 

14 

0.350 

2.20 

0.001745 

E  A?  ■  2.902E-5 

Ei 

a_  -  0.003809 *rad* 

en 

Total  beam  noise  Is: 

en  "  Ke  Z  Aet  sin  (“ie  + 

Where  Ke  ■  1.0 

Beam  noise  is  Injected  into  FLOLS  as  follows: 


Where 

R  -  range 

d  ■  perpendicular  distance  from  nominal  glide  slope 
e  -  angular  deviation  from  nominal  glide  slope 


*od  -  27  ft  at  R  -  7200  ft 


I 
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TABLE  3.  SUMMARY  OF  EXPERIMENTAL  MATRIX 


Conf igu ration  No. 

Student  Code 

Range 

Beam  Noise 

"Loose" 

la,b,c 

0 

Variable 

OFF 

2a,b,c 

0 

Variable 

ON 

3a, b 

0 

Long 

OFF 

4a, b 

0 

Long 

ON 

5a, b 

0 

Medium 

OFF 

6a, b 

0 

Medium 

ON 

7a, b 

0 

Short 

OFF 

8a, b 

0 

Short 

ON 

9a, b 

0 

Variable 

OFF 

10a, b 

0 

Variable 

ON 

"Tight" 

la,b,c 

2 

Variable 

OFF 

2a,b,c 

2 

Variable 

ON 

3a  ,b 

2 

Long 

OFF 

4a, b 

2 

Long 

ON 

5a  ,b 

2 

Medium 

OFF 

6a, b 

2 

Medium 

ON 

7a  ,b 

2 

Short 

OFF 

8a,  b 

2 

Short 

ON 

9a  ,b 

2 

Variable 

OFF 

10a, b 

2 

Variable 

ON 
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that  he  ignored  Che  beam  noise.  On  Che  "StudenC  2"  runs,  Che  piloc  was 
Cold  Co  actively  crack  che  glide  slope,  as  would  be  required  Co  follow  low 
frequency  ship  motions  if  they  were  present. 


ANALYSIS 

The  NIPIP  software  was  used  Co  analyze  a  select  sample  of  Che  runs 
from  Table  3.  The  objectives  of  Che  analysis  were  Co  discern  differences 
in  pilot  control  strategy  due  Co  variation  in  the  following  experimental 
conditions: 


1.  Constant  versus  variable  range 

2.  With  and  without  beam  noise  (surrogate  for  ship  motion) 

3.  Loose  versus  tight  glide  slope  tracking 

The  following  runs  were  selected  to  meet  the  objectives  defined  above: 


Loose 

Tight 

5b. 0 

5b. 2 

6a. 0 

6a. 2 

9b. 0 

9b. 2 

10b. 0 

10b. 2 

Time  histories  for  the  runs  listed  above  are  shown  in  Figs.  2  through  9. 
These  time  histories  were  used  as  inputs  to  Che  NIPIP  software. 

A C  A 0 

The  outputs  of  NIPIP  (i.e.,  Y  .  Y",  and  Y  )  are  summarized  in 

P  P  P 

Figs.  10  through  15.  The  figures  are  arranged  as  follows:  Figs.  10,  11, 


Pitch  Attitude,  9,  deg 
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Altitude,  h,  feet 


T 


Range,  x,  feet 


T 


Axial  Gust,  u  ,  ft/sec 
& 


— i - 1 — 

20  5C 

time,  seconds 


Figure  2  (Concluded) 


10  20  30  4C 

time,  seconds 

Figure  3.  Tine  History  of  Configuration  No.  6a. 0 
(Studen*  0,  Constant  Range,  Ream  Noise  On) 
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0.5 


Glide  Slope  Deviation,  e,  deg 


NAVTRAE 
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Beam  Noise  Configurations 
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Figure  Time  History  of  Clide  Slope  Loop  Describing  Function,  Y  (Jw), 

No  Besa  Noise  Conf l gu rat  ions  ^ 


nf Igurat Ions 


Tlroe  History  of  Crossfeed  Describing  Function,  Y  (jw), 
With  Beam  Noise  Conf lgurat iona  ^ 


History  of  Glide  Slope  Loop  Describing  Function,  Yp(J<*>), 
With  Beam  Noise  Configurations 
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*  0  *  g 

and  12  are  Yp,  Y^,  and  Yp,  respectively,  for  all  of  the  no-beam-noise  runs 
(5b. 0,  9b. 0,  5b. 2  and  9b. 2);  Figs.  13,  14,  and  15  are  the  equivalent  de¬ 
scribing  functions  for  the  runs  with  beam  noise  (6a. 0,  10b. 0,  6a. 2, 
10b. 2) . 

The  time  history  of  each  describing  function  is  presented  as  an  am¬ 
plitude  ratio  (dB)  and  phase  (deg)  at  a  single  frequency,  where 


IvH,.  ■  20  lo*io  lyH 

do 

Thus  the  amplitude  plots  of  Yp(Jw)  are  logarlthalc  scales  (e.g.  ,  a  20  dB 
change  in  |Yp(j<i»)|  is  a  factor  of  10,  a  6  dB  change  is  a  factor  of  2). 
The  phase  plots  Indicate  the  relative  amounts  of  lead  or  lag  between  the 
control  and  the  variable  to  be  controlled. 

While  the  describing  function  at  each  frequency  exists  at  each  com¬ 
puted  time,  its  value  near  the  ciossover-f requency  region  of  each  loop  is 
of  key  importance.  Therefore  we  have  chosen  to  illustrate  Yp  at 
2.0  volts/sec  (about  1/3  Hz)  for  inner  loops  (9)  and  0.5  volts/sec  (about 
0.1  Hz)  for  outer  loops  (x  and  e). 

The  legend  used  to  plot  the  data  In  Figs.  10  through  15  is  as 
follows: 

1.  Open  symbols  were  "loose”  tracking  ("Student  0,"  pilot 
not  following  beam  noise) 

2.  Closed  symbols  are  "tight"  tracking  ("Student  2,"  pilot 
told  to  follow  beam  noise) 

3.  The  symbols  used  were: 

Medivaa  range,  en  ■  0,  Configuration  5 
Medium  range,  en  +  0,  Configuration  6 
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Variable  range,  en  -  0,  Configuration  9 
Variable  range,  en  ^  0,  Configuration  10 


A  sliding  time  window  of  30  seconds  was  used  for  all  of  the  data 
presented  In  Figs.  10  through  15.  Also,  the  first  10  seconds  of  each  run 
was  ignored;  and  NIPIF  waits  5  seconds  before  outputting  the  estimates  of 
Yp(jw).  Thus  there  is  no  output  from  NXPIP  until  t  ■  15  seconds  in  all  of 
the  figures. 

Conclusions  and  observations  based  on  the  time  histories  shown  in 
Figs.  10  through  15  are  presented  in  the  next  subsection. 


COMCLOSIOMS  AMD  OBSERVATIONS 

The  following  conclusions  and  observations  are  based  on  the  somewhat 
small  data  base  presented  in  the  previous  section.  In  order  to  substan¬ 
tiate  these  findings,  we  recommend  that  a  larger  amount  of  data  be 
analzyed.  It  is  not  the  purpose  of  this  report  to  make  broad  statements 
on  pilot  control  techniques;  rather,  the  main  purpose  of  this  report  is  to 
demonstrate  the  kinds  of  conclusions  that  can  be  made  by  using  the  NIPIP 
software. 


Pitch  Loop 

There  were  no  consistent  changes  in  the  pitch  attitude  describing 

A 

function,  Yp(ju),  with  changes  in  experimental  conditions.  Except  for 

Configuration  5b. 2  in  Fig.  10,  the  pilot  tended  to  use  a  slight  amount  of 

lead  (<  Y®(2.0J)  *  0  to  +25  deg),  and  a  gain  between  -20  dB  and  -25  dB. 

*  0 

This  lack  of  trends  in  the  Yp(jw)  data  indicates  that  the  inner-loop  clo¬ 
sure  on  6  ♦  6e  is  relatively  Independent  of  simulated  beam  noise  or 
tightness  of  the  glide  slope  loop.  This  makes  sense  because  the  pilot 
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■use  always  closely  stabilise  the  pitch  attitude  of  the  aircraft  against 
disturbances  by  using  the  elevator,  and  the  tightness  of  this  loop  Is  set 
mainly  by  aircraft  constraints  (Ref.  5). 


Crossfeed 

There  were  some  very  Interesting  and  consistent  trends  In  the  cross- 
feed  Y^(ju>).  For  the  runs  with  beaa  noise  (Fig.  13)  the  pilot  used  a  very 
high  crossfeed  gain  for  the  "Student  2"  runs,  but  virtually  no  crossfeed 
for  the  "Student  Cf "  runs,  where  beaa  noise  was  ignored.  Also  hote  that 
the  gain  was  practically  the  saae  for  both  the  constant  and  the  variable 
range  runs  for  "Student  2"  but  not  for  "Student  0.”  This  Is  consistent 
with  the  experimental  conditions  because,  for  the  "Student  2"  runs,  the 
pilot  was  told  to  actively  track  the  glide  slope;  but,  for  the  "Student  0" 
runs,  the  pilot  admitted  that  he  was  ignoring  the  glide  slope  distur¬ 
bances.  More  about  this  result  will  be  stated  shortly. 

There  was  much  more  variability  and  a  lack  of  consistent  trends 
In  Y*(ju)  data  for  the  no  beam  noise  runs  (Fig.  11).  These  results  indi- 
cate  a  lack  of  consistent  piloting  technique  when  there  is  an  Insufficient 
disturbance  in  glide  slope  deviation.  This  result  has  some  Important 
implications  with  respect  to  pilot  training.  First,  the  pilot-aircraft 
system  must  be  disturbed  from  the  nominal  glide  slope  In  order  for  the 
pilot  to  learn  (and  hence  to  demonstrate)  the  proper  control  techniques. 
The  simulated  beam  noise  used  in  this  experiment  provided  the  required 
disturbance,  but  the  pilot  thought  it  was  unrealistic.  Indeed,  the  pilot 
tried  to  ignore  the  beau  noise  during  the  initial  set  of  runs  ("Student 
0*').  A  more  realistic  glide  slope  disturbance  could  be  easily  provided  by 
using  wind  shears  starting  during  different  portions  of  the  final 
approach.  This  would  not  only  force  the  pilot  to  use  the  proper  control 
technique  but  would  also  acquaint  him  with  the  limitations  of  the  aircraft 
(e.g.,  when  he  should  wave-off  versus  reacquiring  the  glide  slope  and 
continuing  to  a  landing). 
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Gild*  Slope 

There  were  elso  some  Interesting  trends  in  the  flight  path  control 
describing  function,  Y*(Ju»)  .  The  "Student  2"  gains,  |Y^(jh>)|,  were 

consistently  higher  then  the  "Student  0"  gains  (shaded  versus  open  symbols 
in  Figs.  12  and  15),  and  the  variability  in  the  data  with  beam  noise 

(Fig.  15)  was  generally  lower  than  in  the  data  without  beam  noise 

(Fig.  12),  especially  in  the  phase  angles,  <  Y^(ju).  The  higher  gains 
indicate  higher  pilot  workload,  and  the  low  variability  in  the  phase 

angles  Indicate  a  more  consistent  pilot  control  strategy.  These  trends 
are  similar  to  the  Y*  data  discussed  above  and  are  consistent  with  the 
experimental  conditions. 


swmah 

The  NIPIP  software  was  successfully  used  to  analyze  selected  runs 
from  an  aircraft  simulation  of  the  T-2C  on  final  approach  to  an  aircraft 
carrier.  NIPIP  quantifies  the  pilot's  control  strategy  in  terms  of  fre¬ 
quency-domain  describing  functions.  The  parameters  of  these  describing 
functions  are  gains  and  phase  angles,  which  can  be  used  to  infer  pilot 
workload  and  performance.  For  example,  a  high  gain  indicates  tight 
closed-loop  control  (high  bandwidth)  but  also  demands  high  pilot  work¬ 
load.  A  positive  phase  angle  indicates  that  the  pilot  is  generating  lead 

in  his  control  technique  -  a  process  which  also  Increases  workload.  A 

negative  phase  angle  indicates  a  lag  smoothing  control  technique. 

The  NIPIP  results  presented  herein  were  consistent  with  the  various 
experimental  conditions.  The  most  interesting  results  were  the  changes  in 
the  pilot's  describing  functions  with  simulated  glide  slope  disturbances 
(Injected  beam  noise)  and  the  "tight"  versus  "loose"  tracking  runs.  For 
the  "loose"  tracking  runs,  there  was  a  very  low  glide  slope  gain  and  vir¬ 
tually  no  crossfeed  gain.  For  the  "tight"  tracking  runs,  the  pilot 
exhibited  high  glide  slope  and  crossfeed  gains  with  relatively  low 
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variability  in  the  data,  especially  for  the  runs  with  beam  noise.  The 
i4rlication  is  that  adequate  glide  slope  disturbances  must  be  present  in 
order  for  the  pilot  to  demonstrate  his  ability  to  control  the  aircraft 
properly,  as  discussed  above. 

The  main  objective  of  this  experiment  was  fulfilled  in  that  N1PIP  was 
able  to  discern  differences  in  pilot  control  strategy  under  simulated 
flight  conditions.  However,  it  is  recommended  that  further  testing  be 
done  using  a  ntmber  of  pilots  (both  skilled  and  trainees)  under  controlled 
experimental  conditions. 
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